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BED UCTION OF THE PROBLEM OF THREE BODIES. 



BY G. W. HILL, NYACK TUENPIKB, N. Y, 

The object of this article is to find the three differential equations which 
virtually determine the sides of the triangle formed by the three bodies, 
bringing to ones aid all the known finite integrals of the problem. 

Lagrange was the first to treat this question in his Essai sur le Prohlhme 
des Trois Corps (Oewj^res, Tome VI, p. 227); but the formulae lacking sym- 
metry, his editor, Serret, has, in a note, supplied this and pointed out an 
important error into which Otto Hesse, who had investigated this subject, 
(Journal fur die Mathematik, Band LXXIV) had fallen. 

By adopting an orthogonal substitution, at the outset, for reducing the 
number of coordinates from nine to six, we can prevent the masses from 
entering the equations except through the potential function or its deriv- 
atives. In this way symmetry indeed appears to be lost, but there is so great 
a gain in condensation of the formulae, that we can carry out some of the 
eliminations which the previous writers have been content only to indicate. 

Let $, y], Z) f'j v'y C; ^"> 'O" 1 C" be the rectangular coordinates of the 
masses m, m,' m", the expression for the living force will be 

idff ^ 2dp ^ 2de ' 

and A, A', A" being given by the equations 

A' = (6' - $"y + {yf - riJ + (C - Z")\ 
A'^ = (f _ f y + {f -r,y + (C"- c )\ 

A"»={^ -e'f + (=7 -ir + {z -<)\ 

the potential function 

Q vdm . mm' , imm 

ii -r- -\- 



A 'A' ' A" ' 
Without lessening the generality, the origin of coordinates can be put at 
the centre of gravity, when the principle of the conservation of this centre 
will furnish the equations 

mf + mT + m'T = 0, ~| 

my -\- mrj -\- m'rj' = 0, > (1) 

mC + wi'C' + rn'X" = 0. j 
By means of these relations three of the variables can be eliminated and 
the number thus reduced from 9 to 6. This transformation is most ele- 
gantly accomplished by putting 
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^ = ax + ^x', 7) = ay + ^y, C = «2 + /?«', 
f = dx + Px', 7j = dy + Py, C' = «'2 +^z', 
f = a"x + ^'V, :y" = a"y + ^'y', C -= «"^ + /5"^', 

where a, a', a", ^, /9', /3" are six constants which may be so taken that they 

satisfy the five equations 

ma -\- m a =0, ^ 
m/3 + m'i9^ + m^' = 0, 
ma/3 + md^ -\- m'a"^' = 0, 



ma^ -\- md"^ -f m'd'^ = 1, 



(2) 



m/J2 + m:p + m"/3"2 = 1. 

The first two are necessary in order that the equations (1) may be satisfied; 
the third is adopted in order that nothing but squares of differential coeffi- 
cients may occur in the transformed T; and, evidsntly, the last two may be 
adopted without thereby diminishing the generality of the transformation. 
These equations may be solved elegantly in the following manner. Put 
•j/m = A sin ^ cos e, -j/m' = /;; sin ;- sin e, v^m" ^^ h cos y ; 
and adopt the four quantities f, (p, o), co such that 

y'm a = sin f cos {(i)-\-e), |/m ^ = sin f cos {<o'-\-s), 
l/m' d = sin f sin {a)-\-i), j/m' ^ = sin (p' sin (co'+e), 
yw. a = cos ^, ytn p = cos ^ ; 

it is plain that the last two of equations (2) will be satisfied, and the first 
three take the forms 

cos y cos f -\- sin y sin ^ cos w = 0, 

cos y cos (p' + sin ;- sin <p' cos lo' = 0, 

cos ^ cos ^' + sin f sin ^' cos (w — co') = 0. 

Hence if the quadrantal spherical triangle ABC ^ 

is constructed, and the arc AD = y, having any 

arbitrary orientation on the sphere, drawn, and 

BD and CD joined, we shall have (p = BD, f = 

CD, w = ADB, oi = ADC as the general solu- | 

tion of the system of equations (2). 

Then, after substitution. 



rp_ dx^+dy^+dz^ dx^+dy'^'+dz' 



2dt^ 
And if we put 

i7 = /5'-r, 



2df 




f = «" _ a, 






J\ 



v=za? + y* + z\ v' = x'^ + 3/'H«", ^" = xx' + yy'+zz', 
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we have the expressions 

J" =fv + 2fgv" + g'v, 
I* =Pv + yg'v" + g%', 
j"2 =pv+ 2f'g"v" + g"y. 
The equations of motion are now 

df <h' df dy' de~ de' 

^_dQ_ ^^dO^ ^_ dQ 

de dx" de dy" d^~ dz" 

Or, regarding i2 as a function of v, v, v", 

d^x _ ^ dQ , dS d?x' _ 2 , dS dS 

df ^dv^^ dv'" de ^ dv'^y~di/" 

d^z ^ dQ , , dQ d?z' » , dQ , dQ 

= 2z h 2 . • = 2z \- z . 

dp dv ^ dv'" dff dv'^ dv" 

From these, by eliminating the partial derivatives of Q, we obtain 

xd^y — yd^x , x'd^y — y'd^x' f. 

^ "** W ' 

zd^x — xd^z , z'd^x' — x'd^z' ~. 

df "*" W ~~ ' 
yd^z — zd^y , y'd^z — z'd^y' « 

The integrals of which are 

xdy — ydx , x'dy — y'dx , 

dt + dt = ^"^f"' 

zdx — xdz . z'dx' — x'dz , . 

+ J- = KSin/tecosv, 



dt dt 



y dz — zdy . y'dz — z'dy , 



di dt 



sm/tsinv. 



k, ft and v being the arbitrary constants. In addition there is the integral 

of living forces 

r= J2 + A. 

If we put 

„ _ dx'-^-df + (fe» ,_ dx'* + dy* + rfz" „ _ dxdx + dydy + dzdz 
u ^ , u ^ , u ^ , 

it is evident that 
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2 d^ 



dt^ ' 



1 dV _ , _ x'd^x' + y'dY + z'd^z' 

2 df ' df ' 

dH" _.,f_ xd^x' + yd^y' + zdH' 4- ic'd^a; -f- y'd'^y + z'd^j, 



We put moreover 



+ P, 
P- 



xdx' — x'dx -\- ydy' — y'dy -\- zdz' — z'dz 

p __ , 

whence 

xdx' + ydy' -\- zdz' 1 dm" 

dt 2~dt 

x'dx + y 'dy + z'dz 1 dm" 

dt 2"^ 

We have 

d£ __ xd^x' — x'd'^x + yd^y' — y'd^y + zd^z' — z'd'^z 
dt ~ 2de 

By the substitution of the values of -t-^, &c., we obtain 
IdH „ dS , ,,dS 



1 d'^v' , n ,dQ , „dQ 



2 di2 
d(2 



dt)' ' " ciy'" 



— M"=2«"(-^ + 



'(- 



di2\ , , , ,,di2 



dt \dv' dv)^^^ Uv"' 

These equations take simpler forms when the variables are changed as follows, 

yf=l{v + v'), W' = iiv — V'), W" == iv", 

i> = l(u + u'), v' = J(m — u'), v" = u". 
Then they become 



dt^ 
d2w' 

dt^ 
d^w" 

dt^ 



dQ , ,dQ , „ dii 1 n 

.V=W -— + W-^-^ + W T-7-, = — iii, 

aw aw aw 

,dQ , dii 

dw aw 

, ,. „dS , dQ 

^ dw dw" 

do , dii It dii 

dt dvf" dw' 



..(3) 
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If we add to the first of these the equation of living forces v=S-{-h, we get 

-5 = *^ + ^' (^) 

an equation involving only the variables w, w' and w". 

If we square the members of the three equations which constitute the 
principle of conservation of areas, and take their sum, the result will evi- 
dently be an equation which is not changed by a cyclical permutation f the 
letters x, y, z. 

We have the identical equation 
{xdy—ydxy+{zdx — xdzf+iydz—zdyf == {x'+f+z'){dx^+dy^+dz') 

— {xdx+ydy-\-zdzy, 
with the similar equation which is obtained by affixing accents to x, y, z. 

In addition there is the identity 
{xdy — ydx){x'dy' — y'dx')-Y{zdx — xdz){z'dx'-x'dz')-\-{ydz-zdy){y'dz'-z'dy')z= 
{xx' -\-yy' -i-zz')(dxdx' -\-dydy' -\-dzdz')-{xdx' -]-ydy' -\-zdz'){x'dx-{-y'dy-\-z'dz). 

From these equations it will be seen that the equation of the sum of the 
squares takes the form 

vu+v'u'+2v"u" - dv'+do^^+dv"^ ^ 2/,' _ j^.^ 

or, after transforming into terms of the new variables, and for convenience 
writing k for ^k 

p' = 2¥ -\- ■ — -p— ^ 2(wv+wv +w V , .... (5) 

an equation which is symmetrical. 

It is evident now that, since the values of v, v' and v" are known from 
the first three equations of (3), we shall have, as the equations determining 
w, w' and w", (4), (5) and the last of (3), provided we can find a relation 
connecting p with w, w', w", v, u', u" and the differentials of the first three. 

Such a relation can be found in the following manner: — assume the four 
indeterminates X, X', X", X" so that the equations 

= 0, 
= 0, 

= 0, 

are satisfied; and treat the last as if they were equations of condition in the 
method of least squares, that is, multiply the first by x, the second by y, 



xX + cc'X + %X" ^ 


dx -p-,. 
dt^ 


2/X+2/'X+ jz" + 


dy y, 
dt 


dt 


dz -rr/ 

di^ 
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the third by z, and take the sum for a first equation; and so on. In this 
way the normal equations formed from them are 

As the number of these equations exceeds that of those from which they 
were derived, they are not independent, and the determinant, formed from 
the coefficients, vanishes, which is the condition determining p. This equa- 
tion is 

r , v'dv"—v"dv , v'dv"^—%}"dv'dv"^vdv'^~\ 
r 2 , v"dv —vdv" , vdv"^—2v"dvdv"+v'dv^~\ , 
I oF " 2 I I'dv — vdv' , (v'dv-\-vdv')(lv" — v'dvdv — v'dv"'^ | „ „ 
or, expressed in terms of the new variables, 

— 4 [wv — wV — w"v"]/>^ 
I Q Fw'dw" — w"dw' , w"dw — wdw" / , wdw' — w'dw ,r\ 

+T di ^+ — dt — '^—n — ' ]p 

_4 fw <^w'+<^w''+(^w"^ _ dw(/(w^^+w' ') n 
L c^<2 dt dt J 

_. r ^ rfw^ + rfw^"— dw"" _ (^w^d(w'— w" ') ~] , 

■■■■vi r„/- ^w^— <^w^'+<^w"^ Jw^ ^rf(w^— w^ ') n /r _ n 

L df dt dt J ~ "• 

If p is eliminated from this equation by means of its value from (5), we 
shall have an equation of the first degree in p, from which the value of 
this quantity can be derived; this value will be written below. 
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In resum^, we can present our results as follows : — Let the five symbols, 
a, 0,* I)', I)" and p, have the signifiications 

Q = [aw+aV+a"w"]-^+[bw+bV+b"w"]-^+[cw+c'w'+c"w"]-« 

where a, b, c, &c., denote certain functions of ths masses and of a single ar- 
bitrary quantity, 



^% dQ fdQ 

di? (?w dw' 









dV /^_ ^ 
tZi* dw dw" 






"?-^- + ^]^ 



= L 



—2 
—2 
—2 



Xwy'— w'y)^+(w"y— wo")'— (w"y'— w'o'O^ 
dw' -wdw' — w'dw .dw" wdw" — w"dw 

JW" di '"dt M _ 

' dw w'ciw — wdw' .dw" w"dw' — w'dw"~ 

_~dt di ~dt di __ 

~dw w"dw-wdw" I dw' w'dw" — w"dw'~ 



dt 



dt 



dt 



dt 



c.rw"dw' — w'dw" , wdw" — w"dw , , w'oJw — wcfw' ,,1 
\ dt "+ df~~" + dt "J 

Then the differential equntions, which determine w, w' and w", are 



df 



■=^ii + h, 



„2 _ OZ.2 1 dw^ + dw'^ + dw"^ o/ , / / , „ ,n 



dp , dQ 



■w' 



df 

dQ_ 

dw' 



The first and second are of the second order, while the third is of the 
third order. It will be noticed that although the expressions, involved in 
them, are not exactly symmetrical, yet they exhibit some approach to sym- 
metry; and perhaps by a slight change in the variables they may be made 
so. But I have not succeeded in discovering such a transformation. 



*For want of sorts the Greek « is here, and throughout the remainder of this article, writ- 
ten for V. — Compositor. 



